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Shear yielding of glassy polymers is usually described in terms of the pressure-dependent Tresca 
or von Mises yield criteria. We test these criteria against molecular dynamics simulations of defor- 
mation in amorphous polymer glasses under triaxial loading conditions that are difficult to realize in 
experiments. Difficulties and ambiguities in extending several standard definitions of the yield point 
to triaxial loads are described. Two definitions, the maximum and offset octahedral stresses, are 
then used to evaluate the yield stress for a wide range of model parameters. In all cases, the onset 
of shear is consistent with the pressure-modified von Mises criterion, and the pressure coefficient is 
nearly independent of many parameters. Under triaxial tensile loading, the mode of failure changes 
to cavitation. 



I. INTRODUCTION 

The ability to predict the conditions under which a 
material will yield is of great fundamental interest and 
technological importance (jJ-Qj. Over the past three cen- 
turies, a number of yield criteria have been formulated 
that predict whether a combination of stresses on a solid 
will produce irreversible deformation. In this paper, we 
use molecular dynamics simulations to test the applica- 
bility of these criteria to amorphous polymer glasses un- 
der multiaxial loading. 

The yield criterion that is most commonly used 
for polymer glasses is the pressure-modified von Mises 
(pmvM) criterion. The original von Mises criterion jj| 
was based on the assumption that yield occurs in the 
material when the elastic free energy associated with the 
shear deformation -FLhear reaches a critical value. For 
small deformations of an isotropic material, .F s hcar is pro- 
portional to the square of the octahedral or deviatoric 
stress 
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where the Oi denote the three principal stress components 
and p = — (<ti + (72 + f3)/3 is the hydrostatic pressure. If 
one ignores anharmonic effects, the condition of yield at 
constant F s h oa r can be reformulated as yield at a constant 
threshold value of r oct , i.e. 
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For polymers it is often observed p| that tensile pres- 
sure promotes yielding and compressive pressure delays 
yielding, which is not accounted for in Eq. (|^). The pres- 
sure modified von Mises criterion includes a linear 
pressure dependence in the value of the octahedral shear 
stress at yield, 



This pressure modification is motivated by an analogy 
to friction, where the shear stress is linearly related to 
normal pressure, rather than to an energy argument like 
that used to motivate Eq. (||). The coefficient a then 
plays the role of an internal friction coefficient. Eq. (||) is 
also employed to describe the shear yield stress in gran- 
ular materials [||. 

The pressure-modified Tresca criterion (pmT) is also 
often quoted in the context of polymer yielding. Here the 
relevant quantity is assumed to be the maximum shear 
stress 
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The prefactor in Eq. (||) is chosen so that, for the same 
To and a, the pmT and pmvM criteria coincide when any 
two of the three principal stress components are identi- 
cal. In all other cases, the pmT criterion gives a lower 
yield stress. 

Several experimental studies are available that address 
the macroscopic yield behavior of polymers. In these ex- 
periments, the material is typically strained at a con- 
stant strain rate using convenient geometries. The most 
commonly studied stress states are uniaxial tension or 
compression, pure shear, and plane strain compression. 
Several different definitions of the yield point have been 
employed to determine the yield stress. Most authors 
f^PUl(| take the maximum of the stress-strain curve as 
the yield stress. The offset-stress Q is also used, par- 
ticularly when the stress-strain curve exhibits no clear 
maximum. It is given by the intersection of the stress- 
strain curve with a straight line that has the same ini- 
tial slope as the stress-strain curve but is offset on the 
strain-axis by a specified strain, e.g. 0.2% (see Figure 0). 
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The offset stress is motivated by the idea that if the load 
was removed, the sample would relax elastically along 
the straight line to an unrecoverable strain that is equal 
to the offset stress. However, the relaxation is generally 
more complicated. The authors of Ref. [O] argue that 
the true yield point ought to be determined from the 
residual strain measured after unloading. They define 
the yield stress as the smallest stress value that gives a 
nonzero residual strain. 

Raghava et al. JTTJ collected yield data for 
polyvinylchloride (PVC), polystyrene (PS), polymethyl- 
methacrylate (PMMA), and polycarbonate (PC). De- 
spite the fact that different definitions of yield were 
used, all data fit the pmvM criterion |l| with the same 
value of a ~ 0.18. More recently, Quinson et al. JlJ] 
concluded that PMMA is always well described by the 
pmvM criterion, but that PS follows the pmT criterion 
at low temperatures and PC always follows the pmT cri- 
terion. They associated the validity of the pmT criterion 
with the occurrence of shear bands. However, Bubeck 
et al. [[l0| measured a larger number of yield points for 
PC and found that they fit the pmvM criterion. This 
discrepancy may be due to the use of different yield 
point definitions in the two studies. The small number 
of uni- and biaxial stress states considered in experimen- 
tal studies makes it hard to distinguish the two criteria. 
No experimental tests of yield criteria appear to have 
considered triaxial stress states. 

On the simulation side, both Monte Carlo (lj] and 
Molecular Dynamics simulations |l5| ] have been presented 
for uniaxial deformation of polymer glasses. Coarse- 
grained stochastic dynamics of shear flow of local vol- 
umes has also been considered Jig] . These models were 
shown to qualitatively reproduce experimental stress- 
strain curves and their dependence on temperature and 
strain rate. However, we are unaware of any studies that 
address multiaxial stress states and relate them to yield 
criteria. 

This paper is organized as follows. In Section || we 
discuss the details of our model interaction potential for 
polymer glasses and our method for imposing and mea- 
suring stress and strain. Section 111 begins with a com- 
parison of yield stresses calculated from different defini- 
tions of the yield point. We then determine both the 
maximum stress and the offset stress in multiaxial stress 
states and compare the results to the pmvM and pmT 
criteria. In particular, we address the nature of the pres- 
sure dependence of the yield stress and investigate the 
validity of Eqs. (||) and (||) in different physical limits. 



We conclude with a discussion in Section IV 



II. POLYMER MODEL AND SIMULATIONS 

Our model of a polymer glass builds upon extensive 
simulations of polymer melt dynamics p7| . The solid 
consists of linear polymers each containing N beads. 



Pairs of beads of mass m, separated by a distance r, in- 
teract through a truncated Lennard-Jones (LJ) potential 
of the form 
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where uq and d are characteristic energy and length 
scales. Adjacent beads along the chain are coupled by 
the FENE potential @ 
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where Ro = 1.5 d and k = 30uo/d 2 are canonical choices 
that yield realistic melt dynamics \Tj&. Mappings of the 
potential to typical hydrocarbons Il7| give values of uq 
between 25 and 45 meV and d between 0.5 and 1.3 nm. 

The above potentials produce very flexible polymer 
chains, since no third-, or fourth-body potentials are 
present. Melt studies of the bead-spring model have re- 
cently been extended to include the effect of polymer 
rigidity Following these authors, we consider a 

bond-bending potential for each chain, 
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where ri denotes the position of the ith bead along the 
chain, and b characterizes the stiffness. 

The equations of motion are integrated using the 
velocity- Verlet algorithm with a timestep of dt = 
0.0075 tlj, where tlj = y/ md 2 j uq is the characteristic 
time given by the LJ energy and length scales. Peri- 
odic boundary conditions are employed in all directions 
to eliminate edge effects. We consider two system sizes 
of 32768 beads and 262144 beads at two temperatures of 
Th = 0.3uo/fcs and T\ = 0.01 ito/fes. Th is very close to 
the glass transition temperature T g ss 0.35 uo/fce for this 
model pfj| , and T\ represents an effectively athermal sys- 
tem. The temperature is controlled with a Nose-Hoover 
thermostat (thermostat rate 2 t£j ) 

To study the yield behavior, we impose tensile or 
compressive strains e.; on one or more axes of the ini- 
tially cubically symmetric solid at constant strain rates 
of dti/dt = 10~ 4 T^j 1 or less. Figure l] illustrates some of 
the important limiting stress states that are discussed be- 
low. For uni- and biaxial studies of yielding, the stresses 
<Ji in the remaining two or one directions are maintained 
at zero by a Nose-Hoover barostat (barostat rate 0.1 t^j) 
pl| . Note that since tlj ~ 3 ps, the strain rates employed 
are 8-9 orders of magnitude higher than typical experi- 
mental strain rates. Nevertheless by comparing to typi- 
cal sound velocities we can ensure that our strain rates 
are slow enough that stresses can equilibrate across the 
system and thus loading proceeds nearly quasistatically. 
A detailed study of the rate dependence of yielding in 
polymers will be presented elsewhere. 
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FIG. 1. Illustration of some limiting stress states employed 
in this study. For uniaxial simulations 02 and 03 were kept 
at zero using a Nose-Hoover barostat. For biaxial cases 03 
was kept at zero. The controlled strains were increased at a 
constant rate. 

All stresses in this paper are true stresses, i. e. a = 
F/A, where A is the instantaneous cross-sectional area. 

Likewise, we use true or logarithmic strains e = J ; * = 
ln(Z/?o)- In all cases, stresses are averaged over the entire 
simulation cell. 

Our model lets us address a variety of different physi- 
cal situations. The influence of adhesive interactions on 
yielding can be examined by varying the LJ potential 
cutoff distance from r c — 2 1 / 6 d, which gives a purely 
repulsive interaction, to r c — 1.5 d or r c = 2.2 d, which 
include progressively longer ranges of the attractive tail 
of the LJ potential. As a byproduct, the density of the 
solid at zero pressure changes with r c . 

Polymer dynamics in the melt is greatly influenced 
by chain length, particularly when entanglements are 
present. To study the effect of chain length on yield- 
ing, we consider an entangled case N = 256 k 4 x Af e 
and an unentangled case N — 16 ss 1/4 x N e , where 
N e w 60 is the entanglement length of the flexible bead- 
spring model |^2j . The bond-bending potential Vb is also 
varied from the flexible limit (b — Ouq) to a semiflexible 
case (b = 1.5 u ). 



III. RESULTS 

A. Onset of yielding and yield stress definitions 

We begin by considering a uniaxial tension simulation, 
in which the polymer glass is expanded in one direc- 
tion, while zero stress is maintained in the perpendicular 
plane. Figure || shows the resulting stress-strain curve. It 



exhibits the typical features of a polymer glass: linear re- 
sponse is followed by a nonlinear region that bends over 
into a maximum. Note that the maximum is reached 
at at strain of about 6%, which is in very good agree- 
ment with experimental stress strain curves, see e. g. [ ^2[ . 
Once past that maximum, the stress drops slightly and 
the material begins to flow, but these large strains are 
not shown. 
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FIG. 2. Stress-strain curve (dashed line) in uniaxial load- 
ing. The maximum stress is indicated by • and the 0.2% 
offset stress by A. The latter is determined by the intersec- 
tion with the straight line, which has the same initial slope as 
the stress-strain curve, but is displaced by 0.2% on the strain 
axis. Also shown is the residual strain e r as a function of the 
maximum strain during uniaxial (□) or biaxial (o) loading. 
The biaxial results are plotted against the effective strain e e ff 
(Eq. (|)). Here T = T h N = 256, r c = 1.5 d, b = 0u , and 
there are 32768 beads. 

Fig. also illustrates how to obtain the maximum 
and offset yield stresses as described in the introduction. 
Both quantities are easy to identify for uniaxial stress. To 
make contact with the approach of ref. Jl2| , we have also 
determined the amount of residual strain after unloading. 
Squares show the residual strain e r as a function of the 
maximum strain before unloading. Results for e r > 0.2% 
are comparable to the residual strains given in ref. Jl2| . 
These authors identified the yield point with the strain 
where a linear extrapolation of e r vanishes. We have ex- 
tended our studies to much smaller residual strains and 
find no evidence of a threshold strain below which the 
residual strain vanishes. Studies of other glassy systems 
also suggest that at T = a finite unrecoverable strain 
should be present for arbitrarily small initial strains, but 
the theoretical situation at finite T is more complicated 

II- 

Note that in order to obtain 0.2% residual strain, the 
initial extension has to be about 4%. This is about twice 
as large as the yield strain from the 0.2% offset crite- 
rion. The offset criterion assumes elastic recovery along 
the line determined by the elastic moduli of the unper- 
turbed solid. Examination of the recovery curves shows 
that they are non- linear. This confirms that much of the 
nonlinear response under loading to the triangle in Fig. S 
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is due to recoverable, but anharmonic deformations. 

The above discussion raises questions about both the 
applicability of the unrecoverable strain definition of the 
yield point and the motivation for the offset stress def- 
inition. However, since the offset stress is widely used 
in engineering applications, we consider both it and the 
maximum stress definition in the following sections. 



Evaluating the yield stress in multiaxial stress 
states 



Figure ^ shows the principal stress components ai as 
well as the octahedral shear stress r oc t for (a) biaxial and 
(b) triaxial loading. In multiaxial stress states one needs 
to plot the stress (es) against a scalar effective strain e c ff 
constructed from the three principal components ej of the 
strain tensor. This choice is not unique, and influences 
the value obtained for the offset stress. A convenient 
choice that we employ in this paper is 
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where v is Poisson's ratio. This expression is propor- 
tional to the octahedral strain and reduces to e\ under 
uniaxial loading, ft also reduces to the effective strain 
defined in Ref. (T^| for biaxial strains, and allows us to 
treat all strain states on the same footing. 

Fig. |^(a) illustrates the maximum and offset definitions 
of the yield point for biaxial stresses. The maximum 
of the stress-strain curves is straightforward to evaluate 
(solid circles), and both nonzero stress components peak 
at the same strain. By contrast, the offset definition gives 
yield at slightly different strains (solid triangles). This is 
troubling, since a yield point should refer to a unique 
instance in time or strain. Moreover, the offset stress de- 
pends on arbitrary factors in the definition of the effective 
strain and in the amount of offset at the yield point. 

We also evaluated the residual strain for different bi- 
axial loading conditions and added them to Fig. || (dia- 
monds). Interestingly, these residual strains fall onto the 
same curve as the results from uniaxial extension when 
plotted against e c ff . As above, there is no indication of a 
finite threshold below which the residual strain vanishes. 

Fig. ||(b) shows that determining the yield point is even 
more complicated under triaxial loading. For this case, 
only the tensile stress component exhibits a clear maxi- 
mum. The compressive components decrease monotoni- 
cally over the range shown. Moreover, the slope of these 
curves increases in magnitude with strain so that they do 
not intersect the line for the offset criterion. 
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FIG. 3. Nonzero principal stress components in the model 
polymer glass during (a) biaxial and (b) triaxial loading. The 
maximum and offset yield stresses are indicated by • and A. 
respectively. Neither the maximum nor the offset definition 
can be applied to two of the three components in (b) . Also 
shown is the octahedral shear stress r oct (thick dashed line) 
and the yield stresses that result by applying the yield point 
definitions to it. Here T = T\, N = 256, r c = 1.5rf, b = 0, and 
there are 32768 beads. 

In order to make progress for general multiaxial stress 
states, we thus decided to apply both the maximum and 
offset stress definitions not to the three stress components 
individually, but to the octahedral shear stress r oct . This 
is a natural choice given that T oct is the quantity that 
enters the pmvM criterion. The method is illustrated in 
Fig. U and again in Fig. [|, where r oct is plotted against the 
effective strain in a variety of uni-, bi-, and triaxial stress 
states. For the biaxial states considered in experiments 
the yield strains obtained from o\ and <Ji separately are 
nearly the same as that determined from T oct flf|| . Thus 
the value of r^ ct is the same as that constructed from the 
separate yield components. 

In the initial elastic response, stress and strain tensors 
are linearly related by a combination of the elastic mod- 
uli Cij of the solid. The moduli are defined by Cy-e^ = <jj, 
from which we obtain 

a/2 (en + 2c 12 )(cn - c i2 ) 
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by straightforward insertion. Thus a plot of r oct versus 
e e ff should collapse all curves onto a stress-strain curve 
with common initial slope. Fig. [| shows that this condi- 
tion is satisfied for two different temperatures and system 
sizes. All curves start out with the same slope, but they 
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split apart rapidly at higher strains. This reflects the in- 
fluence of pressure on yielding. We will show later that 
this pressure dependence is well described by Eq. (||) for 
both yield point definitions. As in experiments, the yield 
stresses decrease with increasing temperature. 




FIG. 4. Octahedral shear stress r oc t versus effective strain 
in simulations with (a) 32768 and (b) 262144 beads at 
Th = 0.3ito/fc.B and T\ — O.Oluo/fcs for different stress states. 
At low strains, all curves for a given temperature collapse 
onto a line whose slope is determined by the elastic mod- 
uli (Eq. (|io|)). At higher strains the curves split apart, with 
the higher curves corresponding to higher pressures. Straight 
lines are drawn at an offset of 0.2% strain for each tempera- 
ture. Here TV = 256, r c — 1.5 d, and b = Ouo- 

The larger systems in Fig. |](b) exhibit much smoother 
curves, since many more independent yield events con- 
tribute to the average stress response. However, since 
the values of yield stress and strain do not change sig- 
nificantly with size, we use systems with 32768 beads in 
most cases. This speeds up the computations and allows 
a larger parameter space to be explored. Note that each 
curve represents the stress response for one particular 
initial configuration of polymers in the glassy state. We 
have observed variations on the order of 5% in the mag- 
nitude of T y with different 32768 bead configurations, but 
the general features of yielding remained unaltered. 



C. Biaxial yielding 

Experimental tests of yield criteria have generally con- 
sidered uniaxial or biaxial loading conditions. In the bi- 
axial case, the pmvM criterion predicts that the values of 



<7i and 02 at yield should lie on the surface of an ellipse 
in the a\ — U2 - plane. The equation for this ellipse can 
be obtained from Eq. (0) by letting 0-3 = 0: 
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A nonzero value of a breaks the tensile/compressive 
symmetry and shifts the ellipse towards the lower left 
quadrant (where the stress is compressive). The amount 
of the shift increases with a. The pmT criterion predicts 
a yield surface bounded by straight lines. For the same 
To and a, it coincides with the pmvM criterion whenever 
two of the stresses are equal. 
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FIG. 5. Yield points for biaxial stress states. Filled sym- 
bols (•, A) correspond to the yield point at T\ and open sym- 
bols (o, A) to Th. Circles refer to the maximum yield stress 
and triangles to the offset-stress. Data in (a) was obtained 
for 32768 and in (b) for 262144 beads. The ellipses represent 
the pmvM yield surface, Eq. (p"l|), and the dashed lines show 
the pmT yield surface. Statistical error bars are indicated by 
the symbol size. Here N = 256, r c = 1.5 d, and b = Otto. 
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In Fig. we perform an analysis analogous to the ex- 
perimental studies, where we only include biaxial data 
from Fig. [|. We employ reflection symmetry about the 
line <Ti = (72, since the solid is isotropic. We have tested 
this isotropy in selected cases and find variations on the 
order of the symbol sizes. The offset stress is generally 
lower than the maximum stress and becomes very dif- 
ficult to evaluate at Th due to the small magnitude of 
the yield stress and the large thermal noise. The yield 
stresses are larger in the lower left quadrant, indicat- 
ing that yield is pressure-dependent. Also shown are the 
yield surfaces predicted by the two yield criteria. As can 
be seen, the data from both definitions of the yield stress 
are better described by the pmvM criterion (Eq. (fil])) 
than the pmT criterion. The values of To and a can be 
found in Table I (see also the next section). In the fol- 
lowing discussion we focus on the pmvM criterion, which 
provides a better fit for all cases studied. 



D. General stress states 



For general stress states it is simplest to test Eq. (g) 
directly by plotting T% ct versus pressure. This plot lets 
us include triaxial stress states, which could not be dis- 
cussed in Fig. |^. In Fig. ^(a) we show the maximum 
shear yield stresses from all stress states at two tem- 
peratures. In general one would expect Eq. (|J) to hold 
as long as the mode of failure is shear. The linear law 
should break down when cavitation occurs in the poly- 
mer. Dashed lines in Fig. |^(a) indicate the threshold 
for cavitation, which only occurs in our model when all 
three stress components are tensile. For isotropic tensile 
loading r^ ct is rigorously zero and can not be used to 
determine a maximum or offset stress. The endpoint at 
r oct = in Fig. ||(a) was determined from the pressure 
where all three stress components reach a simultaneous 
maximum. The offset stress is not well-defined in this 
limit because e c g also vanishes. 

When yield occurs through shear, the yield stresses 
shown in Fig. |^(a) lie on a straight line as predicted by 
Eq. (^). To test this linear dependence over a wider 
range, we have also performed simulations in which the 
simulation cell was first compressed isotropically to a 
higher pressure. The same range of strains was then im- 
posed again, and the resulting shear stresses were added 
to Fig. g(a). This procedure led to all the points for 
p > 1.5uo/d 3 , which lie on the same straight line as data 
from the uncompressed starting state. Other simulations 
show that this linearity extends to even greater pressures 
of up to 20u /d 3 . 

Fig. H(b) shows the offset yield stresses for two tem- 
peratures. While, the data from the compressed and un- 
compressed starting states each fall on a straight line, it 
is no longer possible to make a common fit through both 
data sets. 
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FIG. 6. Octahedral stress at yield versus pressure at T\ 
(filled symbols) and Th (open symbols). The maximum stress 
is shown in (a) and the offset stress in (b). Solid lines are fits 
to Eq. (|^), and the values for to and a are given in Table I. 
Results for p > 1.5 Mo /d 3 were obtained by applying strains 
to a compressed state (see text). Points where cavitation was 
observed are connected by a dashed line in (a) and are not 
shown in (b). 

Moreover, the values of a tend to be smaller for the 
compressed state than the uncompressed state. From 
Fig. [| one sees that the value of a from the offset crite- 
rion depends on the rate at which stress increases with 
pressure at fixed strain, and the slope of the stress-strain 
curve near the intersection. Both factors change with 
temperature and the choice of offset at yield. For this 
reason, it is not surprising that values of a from the off- 
set criterion show larger variations than those from de- 
termined from the maximum stress. 



E. Effect of adhesive interactions 

Intuitively one would expect that increasing adhesive 
interactions between polymers would lead to a higher 
yield stress. In our model this can be achieved by in- 
creasing r c , which controls how much of the attractive tail 
in the LJ potential is included. We have verified that in- 
creasing r c to 2.2 d produces higher yield stresses and that 
the yield points are consistent with the pmvM criterion. 
Biaxial data lie on an ellipse rather than the straight line 
segments predicted by the pmT criterion (Fig. @). Val- 
ues of To Ct from general stress states fall onto a straight 
line when plotted against pressure (see Fig. below) as 
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long as failure occurs through shear rather than cavita- 
tion. Fit values of To and a are presented in Table I. The 
value of a is the same as for smaller r c within our error- 
bars. The value of tq increases, leading to a higher yield 
stress at all pressures. Increasing r c also suppresses cav- 
itation. The magnitude of the tensile pressure needed to 
produce cavitation under hydrostatic loading increased 
by 2 uo/d 3 at both temperatures. 
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FIG. 7. Maximum shear stress r^ ct minus to at Ti (filled 
symbols) and 31 (open symbols) for three different potential 
cutoffs r c = 1.1225 d (•, o), r c = 1.5 d (A, A) and r c = 2.2 d 
(T,V). To avoid overlap, results for Ti were displaced verti- 
cally by 0.5uo/d 3 . Points where cavitation occurred are not 
included in this plot. Solid lines show best fits to all points 
plotted. 

The opposite limit that we can consider is the purely 
repulsive potential obtained when r c — 2 1 / 6 d. Such a 
solid does not provide any resistance to tensile stresses 
at p = 0, but the same range of strains considered above 
can be applied at higher offset pressures as done in Fig. If]. 
Once again, the results are in good agreement with the 
pmvM criterion and values of To and a are quoted in Ta- 
ble I. The value of To is greatly reduced from those for 
larger r c , but a remains unchanged within our errorbars. 

To illustrate that a is constant, we summarize data for 
the maximum shear yield stress as a function of pressure 
for all three different ranges of the potential in Fig. ^. For 
each case, the value of To is subtracted, and results for 
T\ are displaced vertically by 0.5 uo/d 3 in order to avoid 
overlap with the high temperature data. For each tem- 
perature, all of the data points collapse onto a common 
straight line whose slope is consistent with the indepen- 
dently determined values of a. Thus while To decreases 
rapidly as adhesion is reduced, a remains unchanged. 
Similar results have been obtained in recent studies |2q| 



of static friction as summarized in Sec. IV 



F. Semiflexible and short polymers 

In the previous section we showed that the interaction 
range has a large effect on tq, but does not significantly 



alter a. Other features of the potential have little effect 
on either quantity. 

We first consider semiflexible polymers with a value 
of b = 1.5 m and a cutoff r c — 1.5 d in Eq. (||). This 
makes the polymer conformation more rigid, and one 
would therefore expect increased resistance to deforma- 
tion. Fig. H (a) shows that the pmvM criterion also de- 
scribes yield of this semiflexible polymer. Values of To 
and a from linear fits in Fig. ||(b) tend to be slightly 
higher than the values for flexible chains (Table I), but 
the changes are comparable to statistical uncertainties. 
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FIG. 8. Yield points for (a) biaxial stress states and (b) 
general stress states from simulations with semiflexible chains 
(b = 1.5 uo). As before, filled (open) symbols correspond to 
Ti (Th) and circles (triangles) indicate the maximum (offset) 
yield stress. Solid lines show fits to the pmvM criterion us- 
ing values of a and to given in Table I. Dashed lines con- 
nect points where cavitation occurred. Here N = 256 and 
r c = 1.5 d and there are 32768 beads. 

Finally, we have considered shorter chains of length 
N = 16 to contrast the behavior of unentangled and en- 
tangled (N = 256) chains. The long-range topological 
properties that are crucial for the dynamics of polymer 
melts turn out to be irrelevant for the small strains re- 
quired to initiate yield (^4). Values for To and a do not 
deviate substantially from the results for N = 256, and 
we merely quote them in Table I. 
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IV. SUMMARY AND DISCUSSION 

We have examined the yield stress of amorphous 
glassy polymers under multiaxial loading conditions us- 
ing molecular dynamics simulations. Measurements of 
the residual strain after unloading indicate that imposing 
any finite strain produces some permanent plastic defor- 
mation. Thus no threshold for the onset of plastic yield 
could be identified fll^| . Two other common definitions 
of the yield point, the maximum and offset stress, could 
not be applied separately to the three principal stresses. 
Instead they were applied to the octahedral shear stress 
T oct for a wide range of loading conditions and model 
parameters. 

In all cases studied, and for both yield definitions, 
the pressure-modified von Mises criterion provides an 
adequate description of the shear yield stress. Tradi- 
tional plots of biaxial yield points are well fit by ellipses 
(Eq. (pT|)) and multi-axial results for T^ ct vary linearly 
with pressure. The pressure-modified Tresca criterion 
always shows worse agreement. Although the difference 
may be small from an engineering perspective (< 21%), it 
has implications about the nature of the molecular mech- 
anisms that lead to yield. 

More experimental studies seem to be consistent with 
the pmvM criterion than the pmT criterion, although 
the number of stress states that has been accessed is rel- 
atively small and limited to biaxial cases. One study ]T^ ] 
found that the pmT criterion worked better for poly- 
mers that formed shear bands. Shear banding can be 
influenced by boundary conditions and the method of 
imposing shear. In particular, our periodic boundary 
conditions inhibit shear banding and this may be why 
the pmvM criterion always provided a better fit to the 
yield stress. It is also possible that the failure criterion 
is dependent on the specific form of the atomic interac- 
tions, and that the simple potentials considered here do 
not lead to pmT behavior. These will be fruitful subjects 
for future research. 

The dependence of To and a on potential parameters 
and temperature was explored. The value of a is nearly 
independent of potential parameters, but seems to in- 
crease slightly with increasing temperature. The value 
from the offset stress is larger and more variable than 
the value from the maximum stress. The value of To 
is also relatively insensitive to the chain length N and 
degree of rigidity, indicating that yield is dominated by 
local structure. However, To increases rapidly with in- 
creases in the attractive interactions between molecules 
that determine the cohesive energy. The value of To also 
decreases with increasing temperature. Temperature and 
rate dependence will be discussed in subsequent work. 

Our calculated values of the dimensionless parameter 
a compare well with typical experimental values. Quin- 
son et al. report values for a between 0.14-0.25 for 
PMMA, 0.03-0.12 for PC and 0.18-0.39 for PS. As in 
our simulations, there was a tendency for the value to 



increase with temperature. Bubeck et al. |h]] found val- 
ues between 0.04-0.07 for PC. Rhagava et al. fit all data 
for the above polymers and PVC to a = 0.18, and best 
fit values for individual polymers ranged from 0.13-0.20. 
They used the offset stress to determine a and our num- 
bers for this quantity are in the same range: a — 0.11- 
0.16. Comparing values of To is more complicated be- 
cause of uncertainty in mapping our potential parame- 
ters to specific polymers (see Section ||). Using values 
of uq and d from Ref. (lj] our values of To correspond to 
1 to 50 MPa. This coincides well with the experimental 
range [||jlj] of 5-100 MPa, but quantitative comparisons 
for specific polymers are not possible without more de- 
tailed models of atomic interactions. 

It is interesting to note a connection between the re- 
sults reported here and recent studies of the molecular 
origins of static friction |25|j2^ ]. Macroscopic measure- 
ments of the force needed to initiate sliding can be un- 
derstood as arising from a local yield stress for interfacial 
shear that rises linearly with normal pressure exactly as 
in the pmvM criterion (Eq. [|). Moreover, the results pre- 
sented here parallel results from simulations of the static 
friction due to thin layers of short chain molecules that 
are present on any surface exposed to air |^6| . These au- 
thors find To represents an increase in the effective normal 
pressure due to adhesive interactions. As in our study, To 
increases with the range of interactions, r Cl and the value 
of a is insensitive to many details of the interactions. A 
simple geometrical model for a was developed based on 
the idea that surfaces must lift up over each other to al- 
low interfacial sliding. It will be interesting to determine 
whether a similar geometric model can be developed for 
yield of bulk polymers. 
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TABLE I. Values for to and a obtained from the linear fits 
in Figures ^, ^ as well as the simulations with r c = 2.2 d. Un- 
less otherwise noted, N = 256 and b — 0. Typical statistical 
errors in fits for the maximum yield stress are ±0.02 uo/d 3 for 
ro and ±0.01 for a. Errors are about three times larger for to 
in the purely repulsive case, and for all offset stress results. 



T\ max 
To/a 



Ti offset 
To /a 



Th max 
To/a 



T h offset 
To/a 



Plathe, e-print cond-mat/0005192 
C. Bennemann, W. Paul, K. Binder and B. Dunweg, 



r c = 2 1/6 d 0.06/0.08 -0.42/0.09 

r c = 1.5d 0.72/0.08 0.37/0.15 0.23/0.10 0.11/0.16 

• = 1.5d, comp. 0.72/0.08 0.24/0.11 0.23/0.10 0.15/0.06 

r c = 2.2d 0.83/0.09 0.53/0.16 0.45/0.10 0.28/0.09 

,= 1.5 d, b= 1.5 0.76/0.10 0.44/0.13 0.31/0.15 0.19/0.11 

= 1.5 d, N = 16 0.68/0.07 0.36/0.14 0.20/0.12 0.07/0.06 







